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Heat Transfer Water Spouted 


GHOSH? and OSBERG* 


Heat transfer from column wall water 
spouted sand beds was studied in. diam. steel 
column. The heat transfer rate water flowing 
through spouted bed was lower than when the 
column was empty. The presence bed suppressed 
convection currents, causing large temperature 
gradient from the wall the centre the bed, and 
consequently lowered heat transfer rate. the 
same water flow rate, the spouted bed gave lower 
rate than fixed bed, because large fraction 
the water bypassed the wall region and flowed 
the spout. 


previous communication from this laboratory 
heat transfer from column wall bed spouted with 
air, Klassen and Gishler reported sharp increase 
the rate heat transfer the onset spouting. They 


also obtained some preliminary data wall bed heat Bulk 
transfer water spouted system. Though there was Material Mesh size density density 
indication that the heat transfer rate decreased the 
onset spouting, the data were not sufficient for publica- 94.5 
tion. The investigation was continued and the results are 
reported the present communication. Medium 160 91.5 
Experimental 
The heat transfer column, shown Figure con- 
sisted two concentric steel cylinders: the inner one was 
in. diameter, 3/16 in. thickness, and fitted with 
60° cone the bottom, while the outer one was in. SUPERHEATER 
diameter, and 1/8 in. thickness. The space between the WATER OUTLET 
cylinders was divided two horizontal rings into three 
independent heating compartments. Steam superheated THERMOCOUPLES 
approximately was used for heating. The enthalpy 
the condensate from each compartment was determined 
The column and all pipings etc. were 
heavily insulated make radiation and other heat 
loss corrections unnecessary. 
Water was introduced the bottom the column 
through half inch orifice. The water leaving the column 
oughly mixed tank, and its final temperature was 
determined before discharging into the drain. The 
inlet and outlet temperatures water were determined THERMOCOUPLES 
mercury-in-glass thermometers. 
Eight chromel-alumel thermocouples were embedded 
the walls the inner cylinder the heat transfer CONDENSATE 
apparatus, and five more were inserted into the column PLATE PIPES 
floating thermocouples for exploring the radial temper- MIXER 
ature distribution different levels. 
TRAPS 
+ 
received July 1959. WATER 
2N.R.C. Post-doctorate Fellow. Present address: Department Chemical FROM CALORIMETER 
Engineering, Jadavpur University, Catcutta-32, India. FLOWMETER 
8Division of Applied Chemistry, National Research Council, Ottawa, Ont. - 
Contribution from the Division of Applied Chemistry, National Research _ 
Council, Ottawa, Ont. Issued N.R.C. No. 5394. Figure 1—Heat transfer apparatus. 
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Coarse, medium, and fine sands were used bed 
materials. Physical properties these are given Table 
the beginning series runs the column was 
filled with water and sand poured into from the top 
fill predetermined stationary height. Water flow 
rate and the steam superheater were then adjusted, and 
when steady state was reached the heat transfer rate 
each section was determined calorimetric measure- 
ments, from knowledge the enthalpies steam 
and condensate, and the rate condensation. The total 
heat flow rate for the three sections was checked against 
that independently determined from the rate water 
flow and the initial and final temperatures water. 


TABLE 
PROPERTIES BED MATERIALS 


h 


The bed within the heat transfer column could not 
visually observed. Hence, in. diam. glass column fitted 
with 60° cone the bottom and provided with half 
inch water inlet was used for visual observation the 
behavior the bed different water flow rates. Also, 
plexiglass split column was used observe the shape 
the spout and estimate the relative flow rates water 
through the spout, and through the annulus surrounding 
the spout. 


Data and Discussion 


The heat transfer areas the three sections starting 
from the bottom were respectively 1.70, 1.92, and 1.89 
ft. the bottom section the cone provided 0.84 sq. ft. 
heating surface and the rest the surface was pro- 
vided the cylindrical portion the column. However, 
distinction was made between heat transferred through 
the conical and the cylindrical surfaces while calculating 
the average rates heat transfer per unit area. Various 
bed heights were used, the maximum being inches. 
bed this height covered the bottom and middle sections 
and part the top section the column. The rate 
heat transfer per unit area through the surface covered 
the bed was approximately independent bed height 
under spouting conditions. typical set data will 
discussed the following section. 


The total rate heat transfer the bottom and 
the middle sections taken together, B.t.u./hr., plotted 
against water flow rate Figure These data refer 
bed height inches. The uppermost curve represents 
heat transfer from wall water absence bed, i.e. 
empty column. The presence bed diminished 
the rate heat transfer considerably shown the 
lower curves. The lowermost curve represents heat trans- 
fer bed fine sand, which spouted the lowest 
water flow rate used. The next higher curve represents 
data medium sand which did not start spouting until 
the flow rate represented the point was reached. 
this point the bed remained stationary, and the 
rate heat transferred from the column walls increased 
with increase water flow rate. However, soon 
spouting started, the heat transfer rate decreased. Similiar 
trends were also exhibited coarse sand which spouted 
water flow rate corresponding the point 
Figure may observed that the rate heat trans- 
ferred from column walls was relatively independent 
particle size long the bed was spouting. For static 
bed, however, the heat transfer rate increased with 
creasing particle size. 
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Figure 2—Total heat flow rate water flow rate. 


206 


should noted that the spout created turbulence 
the region immediately above the bed. water flow 
rates were increased, spout height, and induced turbulence 
also increased, that very high flow rates particles 
the top the bed fluidized. This fluidization caused 
increase 25-50% the rate heat transferred from 
the section contact with the top part the bed. For 
the data shown Figure however, the turbulent region 
affected only the upper section, and did not contribute 
the heat transferred the bottom and middle sections 
which are reported. 

Heat transfer data are usually reported terms heat 
transfer coefficients calculated from the equation 


The main difficulty applying this equation the 
spouted bed data the selection At, i.e. the average 
temperature difference between steam and the bed. 
the absence bed convection currents obliterated tem- 
perature gradients except near the walls making com- 
paratively easy determine At. the presence 
bed, however, convection currents were suppressed, and 
the temperature decreased gradually from the column 
wall the center the bed. This evident from Figure 
which shows typical set temperature readings across 
horizontal section. difficult obtain the average 
temperature the bed any cross-section from inte- 
gration the corresponding temperature distribution 
curve since neither the composition the water-sand 
mixture nor the flow rate the mixture can easily 
measured. Also, continuous record temperature 
any point indicated rapid fluctuations. the central por- 
tion the column temperature fluctuations were about 
0.5°C., while near the walls these were more than 
15°C. (The points plotted Figure represent the 
time average temperatures the corresponding points.) 

Because these circumstances, the overall temper- 
ature difference was obtained follows. The inlet water 
temperature, and hence the inlet, i.e. the bottom 
the section under consideration was obtained direct 
measurement. Since the bed remained confined within 
section the column, the total heat transferred through 
the walls was ultimately carried out the section the 
water flowing through bed. The ‘mixing cup temper- 
ature” the water leaving the section, and hence 
the top the section under consideration, was obtained 
calculation from the inlet water temperature, the 
water flow rate, and the rate heat transferred the 
arithmetic average the At’s the bottom 
and the top the section was used Equation (1) for 
calculating the average overall heat transfer coefficient. 


plot against flow rate presented Figure 
using the same basic data illustrated Figure 
Since the At’s were different for different flow conditions 
the two figures, though similar, cannot superimposed. 

The heat transfer coefficient from steam water flow- 
ing steam-jacketed vertical tube can estimated 
from the natural convection equation recommended 
McAdams Colburn and Hougen Values cal- 
culated from these equations were approximately 
those obtained experimentally. Forced equa- 
tions are not applicable this case since and L/D 
are lower than the specified lower limits for which the 
equations apply The higher heat transfer coefficients 
observed the present apparatus were probably due 
the vigorous circulation water induced the high 
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Figure 3—Temperature distribution. 


velocity jet entering the column. These were relatively 
independent flow rates (Figure 4). 

Heat transfer characteristics water spouted bed 
are similar those stirred slurry. overall coeffi- 
cient 150-175 (°F.) has been reported 
for transfer heat from steam 
slurry crystals jacketed cast iron vessel provided 
with double scraper. For water spouted bed the overall 
coefficients were the order 75-100 indicating per- 
formance similar that poorly stirred slurry. 

Heat transfer studies Richardson and Mitson 
water fluidized beds show that the heat transfer rate was 
more than doubled the presence fluidized bed. 
Similar results were also reported Caldas This in- 
crease the rate heat transfer has been ascribed 
reduction the thickness the laminar layer adjacent 
the walls the column, well the extra turbu- 
lence created the particles. water spouted bed, 
however, the heat transfer rate from wall bed con- 
siderably lower even though there scrubbing action 
the particles against the walls. fluidized bed, there 
are temperature gradients across the bed because 
the rapid mixing within the bed, whereas spouted 
bed large temperature gradients are observed (Figure 
showing that mixing and hence convective heat transfer 
within water spouted bed poor. the lack 
convective currents which largely accounts for low rates 
heat transfer water spouted beds. 

analogy with water spouted beds one would also 
expect suppression convection currents 
spouted beds, and would appear anomalous that 
while water spouted bed shows decrease the rate 
heat transfer through the column walls the onset 
spouting, air spouted bed shows considerable increase. 
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This may, however, explained from the difference 
the physical properties air and water spouted beds. 
the air spouted bed the heat capacity per unit volume 
greatly increased the presence solids. For example, 
sand air mixture containing 50% sand volume has 
heat capacity approximately 700 times that air alone. 
known that this increase heat capacity partially 
contributes the rapid heat transfer beds fluidized 
with air Likewise, owing the increased heat 
capacity, air spouted bed considerable amount 
heat transferred from the solids air the circula- 
tion solids the onset spouting, and this more than 
compensates the effect suppression convection cur- 
rents within the bed. result the overall rate heat 
transfer shows increase. the other hand, water 
spouted bed the heat capacity the bed per unit volume 
not very different from that water. Hence, the 
effect suppression convection currents within the 
bed predominates, and the rate heat transfer from 
wall bed considcrably lowered. For the same water 
flow rate, however, the heat transfer rate from wall 
packed bed higher than that spouted bed, because, 
the latter large fraction the total flow, approxi- 
mately 70-80%, passes through the spout. Consequently, 
the flow rate past the heating surface, and hence the heat 
transfer rate, reduced the onset spouting. view 
the inefficiency wall bed heat transfer water 
spouted bed when compared with wall 
transfer fixed fluidized bed under similar water 
flow conditions, attempt was made develop empiri- 
cal equations for correlating the data. 
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Nomenclature 

average heat transfer coefficient 

heat transferred per unit time 
area heat transfer surface 

overall temperature difference 

average overall temperature difference 

number 
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Influence the Presence Vapor Diffusion 


from Wetted Non-Adiabatic Boundary upon the Sensible 


Heat Transfer between Boundary Wall and Gas 


The evaporation liquid from wetted surface 
into gas stream has consequences which are 
fundamental importance the analysis sensible 
heat transfer the boundary. The presence 
relatively high concentration vapor the region 
adjacent the surface influences the physical pro- 
perties the fluid and the vapor concentration 
gradient affects the heat transfer mechanism. 


Experimental equipment was devised measure 
the sensible heat transfer coefficients the inside 
surface porous ceramic tube from which water 
was being evaporated. Series tests were carried 
out for wall conditions over ranges 
the Reynolds Number from 5,600 16,000, tube 
diameter from 1.360 3.300 inches, bulk air tem- 
perature from 150 220°F., and wall vapor pres- 
sures from 0.35 2.20 p.s.i.a. 


Over the range vapor pressures and evapora- 
tion rates investigated, the sensible heat transfer 
coefficients were not observed depart significantly 
from those predicted generally accepted relation- 
ships for the dry wall case. 


evaporation liquid from wetted surface into 
gas stream has consequences which are funda- 
mental importance the analysis sensible heat transfer 
the boundary region. Individually, heat and mass trans- 
fer have been the subjects extensive experimental 
investigation and analytical work but comparatively little 
attention has been paid the influence mass transfer 
the transfer sensible heat when they occur simul- 
taneously the same field. 


Consideration this combination phenomena act- 
ing together suggests several factors importance. The 
presence relatively high concentration vapor 
the region adjacent the wetted surface influences the 


iManuscript received July 31, 1959. 

2Associate Professor, Department of Engineering Science, Ontario Agri- 
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2Assistant Professor of Mechanical Engineering, University of Toronto, 
Toronto, Ont. 

The work reported here was carried out in the laboratories of the 
Mechanical Engineering Department at the University of Toronto and 
was used Hedlin his doctoral thesis that institution. 


physical properties the fluid and, therefore, affects the 
heat transfer mechanism. Secondly, the diffusion mass 
through the boundary region results sensible fluid flow. 
This phenomenon, which was discovered Stefan and 
studied further Meyer and Thoma “), affects the rate 
sensible heat transfer. Thirdly, sensible heat required 
change the evaporated vapor the gas stream temper- 
ature and thus heat source sink exists this region 
with consequent effect sensible heat transfer the 
boundary wall. Finally, the vapor must accelerated 
stream velocity and the resulting inertial force has been 
shown Ackermann increase the thickness the 
boundary layer and therefore reduce the rate 
sensible heat transfer. 


Several investigators have developed expressions for 
sensible heat transfer which incorporate the mass transfer 
effects. Nusselt derived expression for sensible heat 
transfer between evaporation surface and gas stream 
different temperature, combining the expression for 
the secondary convective stream, caused mass transfer, 
with conventional expression for heat transfer dry 
surface. Spalding’s theoretical work heat and mass 
transfer for fuel droplet combustion yields the expression 


Nu* 


where and are, respectively, the Nusselt numbers 
with and without evaporation. For evaporation without 
combustion 


where the specific heat the vapor, and are 
respectively, the gas temperature and the droplet surface 
temperature, and the latent heat vaporization 
the liquid. Godsave presented theoretical expres- 
sions for heat transfer between flame front and the 
evaporating spherical fuel drop which the front encloses. 
These expressions, together with experimental data, are 
used show that, for burning fuel drops, the Nusselt 
number decreases the mass transfer rate increases. 
formulating expression for heat transfer for the case 
condensation mixed vapors, Colburn and Drew 
included the sensible heat contributed cooling the 
vapor condensing temperature. Experimental heat trans- 
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fer results were correlated Ingebo for adiabatic 
evaporation nine different organic liquids from spheres 
using the Nu, Re, and groups the equation. 


The analogy between heat and momentum transfer 
first proposed Reynolds and since extended 
number workers including Von and 
Martinelli facilitates analysis. 


Molecular and eddy diffusion cause transport fluid 
across turbulent stream and, temperature differential 
exists, transfer sensible heat also results. When the 
Prandtl number equals unity, the equation 


erived from the analogy heat and momentum transfer, 
gives values that agree closcly with experimental heat 
transfer results. Since secondary convective stream 
influences the sensible heat transfer for the 
wetted wall case, this concept particular interest and 
leads method correlating heat transfer results. 


For the case fluid stream flowing along dry 
boundary wall, consider the rate diffusion the flow- 
ing fluid, normal the direction stream motion, (asso- 
ciated with momentum transfer the wall) equiv- 
alent sensible flow having velocity the wall. 
The acceleration quantity fluid pu, the average 
stream velocity accomplished unit time and the 
drag per unit area the wall 


but also 


therefore, 


Following McAdams’ recommendation 


heat transfer, 


and for the friction factor 


0.046, 


one obtains 


The transport heat and momentum the hypo- 
thetical sensible flow therefore equivalent that which 
actually results from molecular and eddy diffusion. 
the particular case heat and mass transfer wetted 
wall, sensible flow gas caused the vapor transfer 
and added the gas interchange that occurs for the 
dry wall case. 


The assumption that 


leads to............ 


Gin ‘ 


the constants Equation (12) are evaluated empi- 
rically, the influence previously mentioned factors 
arising from the mass transfer should accounted for 
with the exception the heat required raise the vapor 
temperature that the fluid stream. 


Equipment 

Experimental equipment was designed and constructed 
permit the evaluation the sensible heat transfer the 
inside surface tube through which turbulent flow 
heated air was passing while non-adiabatic evaporation 
water took place from the surface into the air stream. 

system (Figure consisting fan, stilling tank 
with adjustable waste opening, orifice meter, elec- 
trical resistance heaters, and the necessary piping provided 
air controlled rate and temperature. Preceding its 
passage over the transfer surface, the air flowed along 
approach pipe which functioned reduce entry 
effects acceptable level connection (Figure 
provided for uninterrupted flow from the approach pipe 
into porous ceramic tube. The inner surface this tube 
was kept wet the manner described below and the rates 
sensible heat transfer and water evaporation its inner 
surface were measured. 

Porous ceramic, having capillary potential about 
inches, formed the transfer surface. Tubes this mater- 
ial, with length-diameter ratio approximately 1.6 
were obtained rough form and finished lathe 
inside diameter tolerance 0.001 inch. From four 
six these tube sections were placed end end 
each transfer assembly, the end sections served guards 
minimize radiation and conduction from non-equilibrium 
areas the measuring sections the transfer surface. 

was necessary supply water each tube section 
separately, and course, rate equal the rate 
evaporation from its surface. Two sumps were attached 
externally each tube section; one around the top and 
the other around the bottom. During test, water was 
supplied the top sump rate which was kept some- 
what less than the evaporation rate for the tube section. 
Water from this sump kept the upper portion the tube 
section wet and water carried capillary force from the 
bottom sump supplied the remainder 
Teflon sheet 0.005 inch thickness was used form the 
sumps and separate adjacent tube sections. 

Chromel resistance wire was wound around the tube 
sections that heat could supplied controlled 

rate the inner surface the transfer surface 
temperature and indirectly the evaporation rate. The 
rate heat loss the surroundings from the transfer 
assembly was measured enclosing the ceramic tubes 
poly ethylene tube ith thermopile placed across 
the wall and using this tube heat meter. 

Since was not possible make direct measurement 
the temperature the transfer surface without causing 
undesirable roughness, thermocouples were located 
drilling two holes the tube wall parallel the axis 
the tube and the same axial plane. One thermocouple, 
enclosed small glass tube prevent corrosion, was 
placed each hole. extrapolation temperature- 
position plot, was possible determine the temperature 
the inside surface adjacent each pair thermo- 
couples. 

The completed assembly, consisting the components 
described above, referred the transfer assembly 
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(Figure 2). With this equipment was possible measure 
the rate water evaporation from each tube section 
separately, and measuring directly the heat supplied 
the heaters obtain difference the heat transferred 
from the air the transfer surface. 


When test was made, was necessary operate 
the equipment for one two hours before completely 
stable conditions were obtained the transfer surface. 
After establishment equilibrium the water consumption 
was carefully measured for two four hours and this 
value, together with other measured variables, was used 
compute the sensible heat transfer coefficient. 


Analysis the results required the evaluation the 
phy properties the fluid. was necessary con- 
sider the influence, not only temperature, but also 
the vapor the boundary region the viscosity, density, 
specific heat and thermal conductivi 


The problem extending the validity equations 
include cases where fluid properties are not uniform 
the boundary region has been studied number 
investigators. Nusselt showed that particular equa- 
tion has general validity for gases when logarithmic mean 
values the wall and bulk stream properties were used. 
Other results obtained with large temperature differ- 
ences across the boundary layer, were successfully corre- 
lated when arithmetic mean values the bulk and wall 
temperatures were used. the work reported here, arith- 
metic mean values the temperatures and vapor pressures 
were used the basis for calculating the physical pro- 
perties the fluid stream. 


The viscosity was calculated using the relationship 
Buddenberg and Wilke Thermal conductivities were 
obtained from tables, and specific heats and densities were 
calculated assuming that the influence the vapor 
was directly proportional its concentration. 


Difficulties were experienced from the deposition 
foreign material the transfer surface and this often 
resulted the existence vapor pressure that was less 
than the saturation vapor pressure for water the existing 
wall For this reason, was necessary 
determine the vapor pressure the basis other 
measured variables. The mass transfer coefficient was 
calculated and the corresponding vapor pressure differ- 
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Figure 2—Transfer assembly. 


ence was then evaluated using the known rate mass 
transfer. 


Seventy-six observations the heat transfer 
cient were made the ranges the variables shown 


Table 


The quantity (Equation (8)) was calculated for 
each observation. all cases these values were small com- 
pared ranging from —0.00098 0.00143 (Fig- 
ure 3). There were equal numbers positive and negative 
values; the average the positive values was 0.00028 and 
the average the negative values was 


Attempts were made discover systematic influence 
this quantity the heat transfer factor, 


TABLE 


| 
Variable Range Variable 


Range 


0.0159 0.0176 


1912 2066 


0.17 1.61 
Ib. mass/(hr. )(ft.?) 


20,000 110,000 
ft./hr. 


1.25 7.80 


0.35 2.16 


0.0437 0.0471 
mass/(hr. 


160 220°F. 
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3—Graphical representation experimental results. 


statistical analysis, which employed the method least 
squares, was conducted fit curve, based Equation 
(12), these values. 


The variation the Schmidt number was very small 
and was assumed constant the analysis. 


Inspection the resulting equation indicated that 
did not rationally account for significant part the 
variation The same conclusion was reached with 
respect equations other dimensionless groups ob- 
tained from Equation (11). From this was concluded 
that the experimental error accounted for the observed 
variation. This conclusion supported the closeness 
with which the experimental results for all rates eva- 
poration conform Equation (6). The use all the 
observations gives value 0.0226 for the multiplying 
constant the latter equation place the generally 
accepted value 0.023, and this value does not vary 
systematic manner when the experimental results are 
plotted against evaporation rate. 


For number reasons this investigation was confined 
the air-water combination and result limited 
generality. More general conclusions the effect 
mass transfer sensible heat transfer must await similar 
work with several additional liquids varied character- 
istics. Formulation general relationships based 
experimental evidence would require investigations over 
broader ranges evaporation rates and gas-surface tem- 
perature differentials. 


Conclusions 


The influence water vapor transfer sensible heat 
transfer non-adiabatic wall not significant the 


ranges variation the variables investigated. The 
experimental results obtained these ranges are satis- 
factorily correlated equations that are applicable 
similar dry wall cases. 


This result tends support some extent the 
practice workers psychrometric and drying problems 
ignoring the effects heat and mass transfer upon 
one another when they occur together. means 
expected that the assumption independence 
would retain its validity over broader ranges the 
variables, and most particularly not high mass transfer 
rates. 


Nomenclature 


area, ft.? 

specific heat capacity constant pressure, 

(°F.); for bulk conditions; for film conditions. 

inside tube diameter, ft. 

mechanical diffusion, 

friction factor. 

shear stress, force/ft.? 

mass velocity fluid, tube cross 

section 

evaporation rate, tube wall). 

acceleration due gravity. 

conversion factor )(hr.?). 

heat transfer coefficient, based 

logarithmic mean temperature difference. 

thermal conductivity, for film 
conditions. 

temperature, °F.; for transfer surface; bulk tempera- 
ture air flowing over the transfer surface. 

bulk fluid stream velocity, 

equivalent fluid velocity perpendicular the boundary 

wall; for the dry wall case; for flow resulting from 

mass transfer; for the general case, with without mass 

transfer. 

partial pressure air, ft.? 

partial pressure vapor for film conditions, p.s.i.a. 

absolute viscosity, for film conditions. 

kinematic viscosity, hr. 

density, ft.*; for film conditions. 

Nusselt number, 

Prandtl number 

Reynolds number 

Schmidt number 
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The Carbon Dioxide-Hydrogen 


The phase behavior the carbon dioxide-hydro- 
gen sulphide-methane system, previously reported 
pressures 400, 1000, and 1600 40°F. 
and pressures 1000, 1330 and 1600 
160°F. Correlation the experimental data 
based upon equilibrium ratios calculated for the 
individual components. Equilibrium ratios were 
found function system composition, tem- 
perature, and pressure. Deviation the liquid phase 
from ideal solution behavior, indicated calcu- 
lated activity coefficient ratios, is, general, most 
dependent upon system composition the higher 
pressures and lower temperatures. 


complete representation vapor-liquid equi- 

librium the ternary system carbon dioxide- 
hydrogen has been achieved exten- 
sion the study both higher and lower temperatures 
than that previously investigated. Work has been carried 
out temperature nearer the maximum for which two 
phase equilibrium occurs this three component system; 
and also temperature which below the critical 
temperatures both the less volatile components 
that vapor-liquid equilibrium occurs all three the 
binary systems. The present work permits prediction 
the effect pressure, temperature and system composi- 
tion the distribution the components between the 
vapor and liquid phases and the approach ideality 
the solutions. Experimental techniques and ana- 
lytical procedures discussed Part this series were 
followed closely this work. 


Theory 


Phase behavior ternary system may adequately 
described the construction the series isobaric, 
isothermal, triangular diagrams depicting the compositions 
co-existing vapor and liquid phases. The experimental 
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Sulphide-Methane System 


Part II. Phase Behavior 40°F. and 
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phase composition data readily correlated calculation 
equilibrium ratios; while comparison between the 
behavior the real system and the behavior ideal 
solution may made calculation activity coefficient 
ratios. The concepts the equilibrium ratio 
activity coefficient ratio were discussed briefly Part 
this series. 


The activity coefficient ratio defined the equation: 


while the experimental equilibrium ratio defined the 
equation: 


The general form the isobaric, isothermal phase 
equilibrium diagrams for the ternary system may pre- 
dicted from knowledge the pressure-composition 
diagrams for the binary systems. Since the critical tem- 
peratures hydrogen sulphide and carbon dioxide are 
above 40°F., the temperature used this study 
both these components exist liquids under the test 
conditions. Accordingly, vapor-liquid equilibrium occurs 
all three binary systems—methane and drogen sul- 
phide, carbon dioxide and hydrogen sulphide, methane 
and carbon dioxide. Data the methane- -hydrogen sul- 
system 40°F. have been reported Reamer, Sage 
and Lacey while data the carbon dioxide-hy drogen 
sulphide system between the ice point and the circon- 
dentherm have been reported Bierlein and Kay 
Donnelly and Katz have estigated phase equilibrium 
the methane-carbon dioxide system temperature 
29°F. which below the minimum temperature investi- 
gated this work. 


temperature the system above the 
critical temperature both methane and carbon dioxide. 
Consequently two equilibrium does 
not occur the methane-carbon dioxide system. Reamer, 
Sage and Lacey report data for the methane- hydrogen 
sulphide stem Values for the carbon dioxide- 
hydrogen sulphide system this temperature may 
interpolated from the data Bierlein and Kay (3), 
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Experimental Methods 


The variable volume, visual type equilibrium cell and 
apparatus described Robinson and Bailey 
was employed this study. Modification the equip- 
ment permitted extension the temperature range over 
which phase behavior could investigated. the work 
160°, air was circulated over electrical heating coils and 
then through the constant temperature air bath surround- 
ing the equilibrium cell. temperature 40°F. the 
air was circulated the same manner over series 
finned copper tubes carrying glycol which had 
been cooled standard vapor compression type 
refrigeration system. Temperatures were held constant 


within 0.2 degrees both 40°F. and 160°F. 


EXPERIMENTAL PHASE COMPOSITION DATA 


Total System 


TABLE 


Vapor Phase 


Experimental Results 

The experimental data obtained pressures 400, 
1000 and 1600 abs. 40°F. and pressures 
1000, 1330 and 1600 abs. 160°F. are presented 
Table Figures inclusive represent the correspond- 


ing ternary isobaric diagrams. each case the inter- 
sections the dew point and bubble point loci with the 
axes the equilateral triangle represent respectively the 
vapor and liquid phase compositions the binary systems. 
Close agreement between the experimental results and 
published values for the binary systems was obtained. 
Figure illustrates the phase equilibrium pressure 
between the vapor pressures hydrogen sulphide and 
carbon dioxide 40°F. and, hence, the two phase region 


No. 

Temperature 40°F. Pressure 400 abs. 

0.677 0.0 0.323 0.489 0.0 0.511 0.977 0.0 0.023 
A-7 0.081 0.194 0.477 0.392 0.940 0.027 0.033 
A-2 0.660 0.142 0.198 0.354 0.889 0.070 0.041 
A-3 0.636 0.197 0.167 0.429 0.253 0.318 0.832 0.129 0.039 
0.620 0.259 0.121 0.423 0.352 0.225 0.819 0.148 0.033 
A-5 0.574 0.037 0.404 0.450 0.146 0.702 0.280 0.018 
0.518 0.426 0.056 0.416 0.503 0.03 0.635 0.350 0.015 

0.451 0.549 0.0 0.365 0.0 0.597 0.403 0.0 

Temperature 40°F. Pressure 1000 abs. 

0.594 0.0 0.406 0.28 0.0 0.716 0.878 0.0 0.122 
0.125 0.324 0.291 0.139 0.570 0.756 0.090 0.154 
0.473 0.252 0.275 0.475 0.625 0.221 0.154 
0.320 0.392 0.288 0.229 0.353 0.418 0.412 0.414 0.174 
0.218 0.523 0.259 0.179 0.471 0.350 0.284 0.537 0.179 
0.117 0.627 0.256 0.104 0.567 0.329 0.144 0.680 0.176 

0.0 0.744 0.256 0.0 0.655 0.345 0.0 0.867 0.133 

Temperature 40°F. Pressure 1600 abs. 

0.526 0.0 0.474 0.292 0.0 0.708 0.744 0.0 0.256 
C-2 0.444 0.089 0.467 0.311 0.087 0.602 0.643 0.063 0.294 
C-3 0.437 0.155 0.408 0.340 0.142 0.518 0.488 0.135 0.377 
C-4 0.388 0.437 Critical System Composition for 1600 abs. 

D-1 } 0.908 0.0 0.092 0.844 0.0 0.156 0.968 0.0 0.032 
D-4 0.876 0.054 0.070 0.836 0.061 0.103 0.936 0.031 0.033 
D-5 | 0.861 0.093 0.046 0.822 0.118 0.060 0.911 0.063 0.026 
0.851 0.149 0.0 0.811 0.189 0.0 0.891 0.109 0.0 

F-1 | 0.81 0.048 0.141 0.724 0.066 0.210 0.903 0.023 0.074 
0.732 0.109 0.159 0.704 0.126 0.848 0.064 0.088 
F-3 0.703 0.155 0.142 0.699 0.162 0.139 0.770 0.123 0.107 
0.696 0.203 0.101 0.676 0.100 0.718 0.197 0.085 
0.691 0.218 0.090 unstable sample 

E-1 0.784 0.0 0.216 0.735 0.0 0.0 0.159 
E-3 0.769 0.028 0.737 0.027 0.818 0.018 
E-2 0.755 0.051 0.194 Critical System Composition for 1600 


All compositions are mole fractions. 
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C02 


Figure Phase diagram 40°F. and 400 abs. 


extends from the methane-hydrogen sulphide axis the 
carbon dioxide-hydrogen sulphide axis. Figure the 
system pressure between the vapor pressure carbon 
dioxide and the critical pressure the methane-carbon 
dioxide pair and the two phase region extends from the 
methane-hydrogen sulphide axis the methane-carbon 
dioxide axis. Figure illustrates the case where the system 
pressure above the criticial pressure the methane- 
carbon dioxide system but below the critical pressure 
the two phase region has diminished small area along 
the methane-hydrogen sulphide axis. 


similar shifting the three component two phase 
region occurs with increase system pressure 
temperature 160°F. However, the two phase region 
never extends the methane-carbon dioxide axis 
vapor-liquid equilibrium cannot occur 
system methane-carbon dioxide 160°F. 


o7 os o4 


Figure diagram 40°F. and 1600 abs. 
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CH, 


Figure 2—Phase diagram 40°F. and 1000 abs. 


Equilibrium Ratios 


Calculated equilibrium ratios are presented Table 
and the effects total pressure and system composition 
upon the equilibrium ratio are illustrated graphically for 
temperature 160°F. Figures and The vapor 
and liquid compositions reported this table are based 
smooth curve passing through the experimental 
points. Equilibrium ratios were calculated regular 
intervals arbitrary composition parameter, intro- 
duced permit comparison experimental values taken 
for random system compositions each the pressures 


studied. Values zero and unity represent the binary 
systems methane-hydrogen sulphide and carbon dioxide- 


hydrogen sulphide respectively. The curves 
binary systems were used guides drawing curves 
through the experimental data for the ternary system. 


Equilibrium ratios for methane not occur, any 
pressure, for parameter unity since this corresponds 


Data of Sage and Lacey x 
Ths Work 
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Figure 4—Equilibrium ratios for methane 160°F. 
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Vapor Composition 


TABLE 


Liquid Composition 


Equilibrium Ratio Activity Coefficient Ratio 


(mole fraction) 


Temperature 40°F. 


0.106 0.423 0.944 0.020 


.450 0.330 0.866 0.098 


NK O 


Temperature 40°F. 


.290 0.071 0.639 0.821 


150 0.510 0.340 0.231 


0.0 
0.1 
0.2 
0.4 
0.6 


oo 


Temperature 40°F. 


0.708 0.744 0.0 
331 0.134 0.535 0.547 0.114 


0.0 
0.1 
0.2 


Temperature 160°F. 


0.0 0.844 0.0 0.156 0.968 0.0 

0.1 0.840 0.138 0.961 0.006 
0.2 0.837 0.033 0.130 0.953 0.013 

| | 

0.4 0.830 0.068 0.101 0.935 0.031 
0.6 0.823 0.106 0.916 0.056 
0.8 0.816 0.148 0.036 0.899 
1.0 0.811 0.189 0.0 0.891 0.109 


Temperature 160°F. 


0.1 0.728 0.900 0.012 
0.2 0.724 0.885 0.026 
0.3 0.719 0.083 0.198 0.871 0.042 
0.4 0.709 0.116 0.175 0.839 0.066 
0.5 0.700 0.150 0.800 0.100 
0.6 0.186 0.125 0.761 0.139 
0.7 0.675 0.225 0.100 0.718 0.197 


Temperature 160°F. 


0.0 0.735 0.0 0.265 0.841 0.0 
0.1 0.735 0.027 0.238 0.821 
0.20 0.745 0.049 0.205 0.785 0.044 


system which does not contain any methane. Simi- 
larly, equilibrium ratios for carbon dioxide exist for 
parameter zero. 

The equilibrium ratio must have value unity 
the critical pressure and temperature corresponding 
each particular system This illustrated 
versus pressure curves terminating value unity 
for each value the composition parameter. The critical 
pressure for each system was obtained from plot 
versus critical pressure temperature. 40°F. this 
curve was drawn through the published critical pressure 
for the methane-hydrogen sulphide system, the experi- 
mental ternary critical composition 1600 abs., 
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(mole fraction 


Pressure 400 abs. 


0.034 0.510 14.9 2.08 
13.0 13.0 1.15 10.1 1.89 
0.036 0.498 5.30 11.8 1.18 4.11 1.72 
0.036 0.520 2.25 9.15 1.24 1.82 1.35 
0.033 0.558 1.69 6.97 1.28 1.36 1.07 
0.018 0.596 1.56 6.67 1.30 1.28 0.903 


Pressure 1000 abs. 


0.122 0.324 1.10 1.90 
0.138 0.353 1.73 4.63 1.21 2.30 1.57 
0.150 0.388 1.45 1.91 1.31 
0.179 0.649 0.864 1.90 2.21 0.623 


Pressure 1600 abs. 


0.286 0.460 2.19 1.77 2.02 0.905 
0.339 0.605 1.18 1.58 1.79 0.690 
0.437 1.00 1.00 1.00 3.84 1.49 0.476 


Pressure 1000 abs. 


0.033 0.785 1.05 1.67 0.465 
0.034 0.879 | 2.40 3.82 1.05 1.54 | 0.421 
| 
0.034 0.888 2.97 1.05 0.325 
0.028 0.898 1.90 1.07 0.275 
0.017 0.908 1.74 1.07 1.05 
0.0 | 0.910 1:73 0.00 1.08 1.05 | - 
Pressure 1330 abs. 


0.088 2.77 1.07 1.57 0.39 
0.089 0.817 2.15 2.47 1.08 1.45 0.35 
0.095 0.845 1.75 1.18 0.26 
0.100 0.875 1.50 1.01 0.21 
0.100 0.905 0.91 0.18 
0.085 0.940 1.24 1.18 1.24 0.84 0.17 


Pressure 1600 Ibs abs. 


0.159 0.873 1.67 0.281 
0.160 0.895 1.42 1.49 1.19 0.235 
0.171 1.20 1.35 0.933 0.202 

1.43 0.840 0.168 


0.194 1.00 1.00 1.00 


and estimated critical pressure and composition for the 
methane-carbon dioxide system, based 
data for temperature 160°F., the curve 
was drawn through the published pressure for the 
methane-hy drogen sulphide system, the experimental 
ternary composition 1330 and 1600 abs., 
and the published critical pressure for the carbon dioxide- 
hydrogen sulphide system. 


Activity Coefficient Ratios 
Activity coefficient ratios were calculated regular 
intrevals the composition parameter from smoothed 
experimental equilibrium ratio data. The required theo- 
—, were computed trom the 


retical equilibrium ratios, 


| 
0 
0 


Dore of Berien ond Koy @ 
Tes Work 


1100 1200 1300 
PRESSURE - ibs /in® obs 


Figure ratios for carbon dioxide 160°F. 


properties the pure components using the technique 
discussed Part this study. 

The critical temperature methane below the 
system temperature both 40°F. and 160°F., and hence 


the vapor pressure pure methane the temper- 
ature imaginary. extrapolation known vapor 
pressure data, values 4100 and 9300 abs. were 
obtained for temperatures 40°F. and 160°F. respectively. 
Pressure-volume data for methane 160°F. were obtained 
directly from the literature while those for 40°F. were 


obtained interpolation published values for temper- 
atures 32°F. 70°F. and 100°F. 

Carbon dioxide has critical temperature which 
below 160°F. but above 40°F. the higher temperature 
value 2350 abs. for the vapor pressure was 
obtained extrapolation known values and pressure- 
volume data were obtained interpolation published 
data the lower temperature, pressure-volume data 
for pressures below the vapor pressure were obtained 
interpolation published values 0°F. and 

Both temperatures the investigation were below the 
critical temperature hydrogen sulphide and hence vapor 
pressure data were available. Pressure-volume data, 
the vapor pressure, were obtained directly from the liter- 
ature for temperatures 40°F. and 
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Figure 7—Activity coefficient ratios for hydrogen sulphide 
40°F. 
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Figure 6—Equilibrium ratios for hydrogen sulphide 
160°F. 


Calculated activity coefficient ratios, regular inter- 
vals the composition parameter, are presented Table 
and the effect pressure and system composition 
illustrated for hydrogen sulphide Figure and 
for carbon dioxide 160° Figure The activity 
coefficient versus composition parameter curves are ter- 
minated the maximum value the composition para- 
meter corresponding system which vapor-liquid 
equilibrium occurs the conditions pressure and 
temperature specified. The calculations indicate that for 
all three components the volatility decreased, relative 
that predicted ideal solution laws, increase 
temperature. However, increase pressure decreases 
the volatility the more volatile components 
creases the volatility the less volatile component. 


ACTIVITY COEFFICIENT RATIO 


Figure 8—Activity ratios for carbon dioxide 


160°F. 
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Conclusions 


methane system forms co-existing vapor 
phases over pressure range from 160 abs., the 
vapor pressure hydrogen sulphide, 1949 abs., 
the critical pressure the methane-hydrogen sulphide 
system. Systems having wide range compositions 
form two phases between these two pressure extremes. 
160°F., the pressure range over which vapor-liquid 
equilibrium occurs more limited, extending from 779 
abs. 1660 abs., and the overall system 
must rich drogen sulphide all cases. locus 
ternary critical points exists 40°F. from the binary 
critical point for the methane-carbon dioxide system, 
approximately 1240 abs. the binary critical 
point for the methane-hydrogen sulphide system 1949 
abs. 160°F., the locus ternary ‘critical points 
extends binary critical point for the carbon 
the binary critical point for the methane- -hydrogen 
phide system 1660 abs. 
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Nomenclature 


composition parameter defined the ratio the mole 
fraction carbon dioxide the vapor phase the sum 
the mole fractions methane and carbon dioxide the 
vapor phase, dimensionless. 
fugacity pure component pressure equal its own 
vapor pressure the system temperature, atm. 
fugacity pure component vapor the pressure 
and temperature the atm. 

equilibrium ratio for component dimensionless. 
mole fraction component the liquid dimen 
sionless. 
mole fraction component the vapor phase, dimen 
sionless. 

activity coefficient ratio for component dimensionless. 
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special case liquid-liquid extraction been 
considered where the solute absorption reacts 
form dimer through rapid second order reversible 
reaction with the solvent. 


dimer and monomer have different diffusivities 
(in the same solvent) the rate absorption becomes 
dependent the ratio these diffusivities and the 
isotherm governing the reaction. 


The non linear differential equation governing 
the absorption was solved numerically with the aid 
digital computer for the case resistance 
one phase only. 

Correlations based these solutions were then 
used formulating rate equation. Methods were 
suggested for the determination the quantities 
occurring the “driving force” the predicted rate 
equation. 


liquid-liquid extraction, interphase mass transfer rates 
have been represented generally the rate equation 


where mass transfer coefficient, the inter- 
facial area, and (c* over-all driving force. 


* 


the driving force represents the concentration the 
phase, from which solute transferring, 
equivalent “equilibrium” concentration extract phase 
concentration units. 

one considers transfer between two stirred liquids 
with known interfacial area apparatus such that 
used Gordon and Sherwood the 
transfer rate equal the rate accumulation the 
extract phase 


(la) 


where the volume the extract phase. 
From (1) and apparent that the mass transfer 
coefficient K,, could calculated from the slope 


against time 


graph 
5 } 


~ 
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typical plot such data shown Figure for the 
batchwise extraction propionic acid from water into 
ether. Since the integration accounts for change 
volume during the transfer process, the marked curva- 
ture the plot indicates that perhaps the simple driving 
force (c* used (1) not correct. should 
noted that for the system Figure the distribution 
coefficient was constant over the range concentrations 
the figure. 


well known that when organic acids 
amides are transferred between water and organic sol- 
vent there will molecular associations the organic 
phase and dissociations the aqueous phase. the 
purpose this paper examine how such effects may 
require modification the driving force for the mass 
transfer such solutes interface such sketched 
Figure 


Theoretical Considerations 


the controlling partial differential equation 
Considering solute, HA, diffusing organic phase 

will assumed that there instantaneous equili- 

brium between the monomer and dimer any point: 


The molecules will diffuse according Fick’s law: 

Ac 


@ 
o 
n > 
sore Propionic acid - iso propyl ether- water 


Time (minutes) 


166-126 46 160 25 
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Figure 1—Test driving (Volume varia- 
tion taken into 
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Figure 2—Schematic representation interface. 


where and are the diffusivities and and are the 
point concentrations monomer and dimer respectively. 

The forms Fick’s second law, obtained the normal 
manner, are 


values has been treated extensively Crank and 
others the present authors have not found Equation 
(6) solved the literature. 


Solution the partial differential equation 

shall consider only the simplest transfer situation 
illustrated Figure that is, transfer into essentially 
infinite extent organic phase from interface 
which there constant equilibrium value c,*. the 
initial uniform concentration c,° and one defines 


then Equation (6) becomes 


solved for the boundary conditions 


Utilization the Boltzman transformation (6): 


Here represents fast but finite production dimer 
and should not confused with rate constants that may 
associated with Equation (2). 
Corresponding the semi-instantaneous equilibrium 
described Equation (2) the equilibrium constant 


defined by: 


Using from (4b) and from (5) permits (4a) 
written 


1 0c, 4k, + 0c, 1 Oc, 
where 


should pointed out that the analytical concen- 
tration used 
and define 


then Equation (6) would revert the ordinary form 
Fick’s second law 


and the diffusivities monomer and dimer were equal 
there would difficulty analysing the diffusion 
associated solute. 

Although the solution Equation (8) with variable 
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converts (10) the ordinary differential equation 


and the boundary conditions (11) 


trial solution (13) the Beckman EASE 


analogue computer, initial values and were 


selected (the latter arbitrarily) for the computation. 
Solutions satisfying the second boundary condition (14b) 
could not obtained, however. 

This may have been due the introduction the 
squared relationship into the differential equa- 
tion which provided the possibility having negative 
well positive values which satisfy the differential 
equation and initial conditions. 

numerical solution was found outlined below, 
incorporating the boundary condition 


(since there zero flux zero concentration), 
manner superficially similar that used Wagner 
Corresponding some very small arbitrary value 


there exists value y,. the transformations 

Vo 


4K.ci* ( 
and 
Cy 
x 
Vo" 
219 


(13) becomes 


where 


Equation (17) can simplified even further the 
transformation 


(20) 


For the right side, (20) can integrated 
yield 


give 


having sufficiently small that (15) holds, can 
shown that the integration constant (21) zero. 
follows finally from (21) that 


This latter equation provides initial slope for the 
solution Equation (17) the Runge Kutta equations 
modified this the following steps 
were taken using IBM 650 computer: 

arbitrarily. 

(b) suitable value was assumed. first 

approximation the value obtained for diffusion 
(c) The values 


using (22) and (18). 

(d) The curve against was computed from 


sumption was necessary. 


agreement was found new as- 


Typical computation results 

Figure shows typical curves obtained the above 
procedure for the and values shown. Figure these 
are converted the more useful curves against 
and are compared with the curve for diffusion simple 
molecule diffusivity D,. 

The areas under the curves Figure also provide 
values the total amount solute that has diffused 
any time This can compared with the amount that 
has crossed the interface until this time. 


Thus, total moles absorbed 


c-8. 


C-8.0-995 


0-01 


Figure 3—Numerical solution differential equation for 


From Equation (23) find: 


Equation (24) may used check the internal 
consistency the numerical solutions obtained. 

the preceding derivation the flux the interface 
was written as 


Oc) Os 


whereas probable that only unassociated molecules 
such cross the interface. One could reason that physically 
the equilibrium between dimer and monomer assumed 
this paper exists the immediate neighbourhood 
the interface. the derivation the differential equation 
the isotherm was assumed hold for all 


including and the flux 2p, must therefore 
Ox 


included Equation (23). 
Table are shown typical results obtained from the 
application the area test discussed previously. 


Transfer rates interface 


Let now see how the above computations can 
applied transfer rates interface. 


For the differential Equation (13) may solved 


variables may separated. The solution is: 


dy yuo a 1 


TABLE 


will reduce the equation form which the 


TEST 


3.7 578 581 
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where 


Table 


Using 


R(t) 


where 


and 
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Figure 4—Numerical solution differential equation. 


For other values the following empirical correlation was 
found fit well 


Figure shows plot the above correlation and 
Table shows the correlation for Equation (26a). 


where 


The transfer rate the interface 


Using (26) where this rate expression becomes 


now analytical concentrations used, then (28) can 
written 


Ne 


where 
and 
TABLE 
2.0 0.422 0.444 
3.7 0.042 0.042 


aor+i 


Figure 5—Initial slope correlation. 


For the more general case where the rate 
equation corresponding (29) 


where 
(32b) 
and 


The average rate absorption for short period 
exposure (t,) can now written down using either the 
Higbie Danckwerts method. 


te 


te 4k, \ Tt. 


or 


2 
© 
rey 
7 8 9 10 12 


Figure 6—Distribution plot: propionic 
ether-water (25°C.). 
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06 
0-4 (0) 2.0 


order evaluate the “driving force” (p* 
mine the quantities and with some accuracy. 
Determination The equilibrium constant may 
determined from distribution data simple manner: 
Consider the distribution the acid between water 
(phase and suitable organic liquid (phase 1). 
linear distribution relationship exists between 
acid the two phases (e.g. between simple mole- 
cules, neither associated nor dissociated) this may 
written 


where and are distribution coefficients defined 
this equation. the organic phase 


(cy) (s) 


Analytical concentration 


Hence the distribution relationship becomes 


For weak acids (36) may arranged the form 


= 30a 


This equation slightly different form, was proposed 
Davies al. These authors found their equation 
predict satisfactorily the distribution acetic 
pionic acids between water and number organic 
solvents concentrations IN, 

found from the intercept and slope the plot linear. 
Figure shown the distribution data for the system 
propionic ether-water plotted according 
(36a). 


Association constants obtained this manner should 
accepted with caution, however. Lassettre 
pointed out that hydration the organic phase may lead 
erroneous conclusions regarding equilibrium constants 
this effect not taken into account. 


Determination and has been pointed out pre- 
(7b)), thar the diffusion model 
predicts concentration- dependent diffusion coefficient 
analytical concentration, for the organic phase. 
The “effective” diffusion coefficient was shown 


1) + (37 
Vil + 
the differential diffusion coefficient plotted 
against concentration the model predicts curve the 
form shown Figure 
values and are then given the 
cepr infinite dilution and the steady value high 
concentration respectively. accurate values for can 
determined very concentrations value for 
found noting that 


C=. 


Figure 7—Predicted diffusion curve. 


4K.(D, D>) 


experimental study the diffusion Pamelia 
oil acetone Kaneko the authors reported 
curve corresponding that Figure rough 
computation (by us) from their curve gave value for 
the order 3.7. These authors stated: “The fact that 
there definite relationship between and seems 
indicate that the diffusion associated molecules 
analogous homogeneous diffusion concentration 
dependent system.” (Here concentration gm./100 
cc. solution. 


Summary 

For the absorption compound into organic 
liquid which forms dimer, rate equation has been 
proposed. 

the hypothesis that associated 
species have different but constant diffusivities, these are 
defined Equations (3), diffusion equation was 
derived and numerical solutions were obtained with the 
aid digital computer for the case constant inter- 
facial concentration equal the saturation value, c*. The 
values for initial slope the diffusion curves thus obtained 
were used correlation involving the dimensionless 
constants and (or value the order 0.1 
0.8, and for 1.0, 2.0, 3.7. The results this corre- 
lation were used predict rate equation, using either 
the Higbie Danckwerts mechanism. 

These equations predicted “driving force” un- 
usual complexity which the “bulk” and saturation con- 
centrations and were combined with such quan- 

For the determination these quantities, methods have 
been suggested involving distribution measurements and 
determination differential diffusion coefficient 
function concentration. 


will seen from the final rate equations that the 


deviations from the “normal” driving force, 


further diffusion measurements lend support 
prediction such Equation (37), the quantity may 
unusual importance and felt that the rate 
equations proposed above, although derived for 
limited boundary conditions may become very useful 
starting point for correlations absorption 
tion rates reactions the type discussed occur. 
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Nomenclature 
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(dimensionless 
interface area, cm? 
r 
(dimensionless 
concentration unassociated species, moles HA/litre 
equilibrium concentration 
saturation concentration unassociated species moles 


bulk concentration unassociated species, moles 


(dimensionless 


(dimensionless 
* 

C1 


° 
(dimensionless 

total (analvtical) concentration moles HA/litre 
diffusion coefficient for unassociated species, 
diffusion coefficient for associated species, 
effective differential diffusion coefficient, 
equilibrium constant for association reaction, litres/mole 
mass cm/sec. 

Yo 
a+6 
Yo? 

distribution coefficient 


Danckwerts surface renewal factor 
t = time, secs. 
coordinate direction diffusion (cm 

x 
Boltzmann variable (dimensionless 
=(C 
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Fundamental Aspects Flow 


Part III. Accelerated Motion Particle 


analysis the extensive literature non- 
steady drag forces supports the correlation the 
data means total drag coefficient, which 
appears function the Reynolds Number 
and reduced time parameter which related 
the number particle diameters traversed since the 
initiation the motion. The added mass concept 
shown both completely inadequate 
retically unsound. increase wake turbulence 
resulting from Reynolds Number increases from 
the surface roughness seems dimi- 
nish the acceleration effects. 


Fundamental studies the flow fields around 
blunt bodies reveal the extreme complexity the 
phenomena occurring during acceleration. Explana- 
tions offered for the characteristics the non-steady 
drag coefficient behavior are shown frequently 
variances with these findings. 


flow fields around sphere very low Reynolds 

Numbers which were discussed Part 
series and the wake phenomena described Part 
applied strictly the case aerodynamically-smooth 
spherical particle, moving infinite, turbulence-free 
fluid. Although these idealized conditions are seldom met 
actual conveying systems contacting 
operations, knowledge the fundamental concepts 
underlying the ideal motion essential for the proper 
understanding and evaluation the complicating effects 
commonly encountered chemical practice. 
The rest this series will devoted detailed study 
these complicating factors, concluding with the appli- 
cation this information the study multiparticle 
systems. The present paper concerned with the effects 
acceleration particulate motion. 


has been common practice assume that the stand- 
ard curve for the drag coefficient function the 
Reynolds Number may applied particles undergoing 
accelerated decelerated motion fluid. Lapple and 
Shepherd have ignored any possible unsteady motion 
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effects the development their equations for the 
calculation particle trajectories, and the 
sideration applies the work Gilbert, Davis and 
Altman which was concerned with the relative velo- 
city solids entrained linearly-accelerated combustion 
gases. Khudiakow and Weidner (6) have also used this 
simplification the analysis multiparticle data, 
and acceleration complications have not been considered 
most similar papers dealing with multiparticle solids- 
gas flow. 

Recent experimental studies these laboratories have 
illustrated the inability the standard drag coefficient 
curve predict the motion decelerating droplets 
spray dryer atomized- reactor and 
the motion accelerating particles flash-dryer 
thus preventing the use particulate heat and mass 
transfer information provide generalized design corre- 
lations. 

The discussion which follows deals first with experi- 
mental observations particles undergoing accelerated 
motion together with the superficial methods used 
correlate these data, and then with the limited funda- 
mental information presently available concerning the 
acceleration phenomena. 


ACCELERATED MOTION PARTICLE 
FLUID 


Macroscopic Experimental Evidence 

1826, Bessel reported experiments with spheri- 
cal pendulum moving air and water, and found 
that account for the observed behavior the con- 
ventional equations motion so-called “added mass” 
M,, which was proportional the mass the fluid dis- 
placed the body, would have added that 
the sphere. The constant proportionality was found 
0.9 and 0.6 for the case air and water respectively. 
Others working with spheres oscillating liquids and 
gases also found this concept useful 
explaining the apparent acceleration effects indicated 
their data. Applied the case particle accelerated 
rectilinear motion through still fluid, the equation 
forces the particle would give: 
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where the body force due gravity, the particle 
mass, and the absolute particle acceleration and 
respectiv rely, the fluid density, the particle 
cross-sectional area, the particle volume, and and 
the steady state fluid drag force and drag coefficient 
respectively—these last two terms are assumed 
analysis unaffected the acceleration. The added 
mass concept enjoyed varying degrees popularity but, 
will shown, completely inadequate since 
greatly underestimates the complexity the flow. 


Equation (1) can rew ritten form which would 
make more plausible inference the nature the 
acceleration effect, indicating that acts such manner 
change the fluid resistance from the steady state 
value different one which becomes func- 
tion the acceleration. Furthermore, assumed that 
can related the acceleration and velocity terms 
the expression: 


where and are empirical coefficients, that 
balance forces would give: 


third—and possibly the best—method accounting for 
the altered drag due the acceleration expressing 
terms modified drag coefficient, defined by: 


One the earliest studies rectilinear motion which 
hinted possible acceleration effects was reported 1900 


cmysec? 


cm/sec? 


1000 2000 3000 4000 5000 


Figure coefficients accelerating spheres falling 
(Lunnon). 


Allen who noted that the drag coefficients 
steel spheres slowly accelerating water were increased 
from percent. Later Schmidt recording 
the movement wax spheres water and balloons 
air, found their accelerations slower than would 
predicted from steady flow drag coefficients, his 
results also suggested that the values were function 
the magnitude the acceleration. Schmidt made the 
interesting observation that eddy formation the rear 
spheres undergoing acceleration would occur lower 
Numbers than would normally appear steady 
motion. This would sure increase the form drag 
the body, and would not accounted for the “added 
mass” concept. 


Increased values during acceleration are also 
indicated the work Magnan and St. Lague and 
Cook for spheres moving water, and 

Hesselberg and Birkeland Relf and Jones (23) 
and Hirsch for spheres moving air. Cook, Hirsch, 
and Relf and Jones reported added mass coefficients which 
varied from 0.46 0.83, and seemed affected 
the acceleration, the fluid density and the past history 
the motion but the exact inter- relationships were obscured 
experimental errors. Acceleration drag phenomena are 
very difficult study since continuous time-distance 
record the particle movement required with sufficient 
precision give accurate values the second derivative 
criterion which has rarely been attained the investi- 
gations thus far considered. 

series investigations started 1924 and which 
are noted for their excellence technique, Lunnon 
(25, 26, 27, 28, studied the dynamics spheres var- 
ied density falling air and water; for the former, the 
precision was such that the maximum deviation between 
similar trials the air was the order 0.002 second for 
distances meters thereby furnishing usable accele- 
ration data. Lunnon found his results approximated 
Equation (2), although the initial stages motion 
the resistance was higher and the final stages lower than 
the average resistance throughout the acceleration region. 
Lunnon was able construct smooth curves through his 
results for air when plotted them co-ordinates which 
were directly proportional and Re. The authors 
this survey have replotted some Lunnon’s data more 
conventional manner, using and the ordinate 
and abscissa respectively (Figure 1), and these show that 
even moderate accelerations can have marked effects 
the drag coefficient, depending the particle Reynolds 
Number regime. These experimental results may not 
too general, since they were all obtained from trials 
which the spheres were equal Lunnon’s data 
for water were not accurate presented 
such analysis, but they indicate that will 
affected rably the density. 

The “added mass” analysis when applied Lunnon’s 
water trials show variations the coefficient from 
2.0, which are not excessively removed from the 
values quoted previous authors. However, the data for 
give values which are high 1,200—results which 
the added mass concept utterly fails explain. 

Lunnon feels that the increase resistance due 
acceleration may result from the rapid change stream 
lines which could favorable the production eddy 
motion. points out that this the case, negative 
acceleration should also cause increase that 
Equation (2) must modified, since predicts decrease 
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resistance with deceleration. Cp, indeed reflects the 
buildup decay vortex patterns, then the rate 
change the acceleration must taken into considera- 
tion. The accuracy required for this type investigation 
has not hitherto been attained, but has now been “made 
possible radioactive tracer method combined with 
digital computer analy sis (31), 

1937, Vagner attempted extend Lunnon’s 
data orking with series liquids cov ering wide 
density spectrum. Although much care was taken 
develop suitable photographic technique, could not 
measure the drag coefficient directly, but had fall back 
the added mass type relationship given Equa- 
tion (1). superimposed his time-distance data 
curves drawn assuming various values, and 
this way found coefficients varying from 0.5 1.5. Al- 
though results show inverse relationship between 
diameter, the gravitational constant, and the fluid 
viscosity, the data not establish this quantitative 
manner. 

Following the work Vagner, the study the 
acceleration effect the drag lay more 
results obtained with droplets air Laws 
and correlated which again suggested 
increase the drag coefficient. The post-war investiga- 
tions the acceleration effect can divided into three 
very general categories, depending whether the bodies 
more relative the fluid—(a) rectilinear but variable 
acceleration; (b) rectilinear and constant acceleration; 
and (c) with oscillatory motion. 


Rectilinear, Variable Acceleration 


Iverson and Balent carefully studied the accelera- 
tion circular discs towed fine eighted cable and 
set perpendicular the motion, and the added mass con- 
cept completely failed explain the large dependence 
the fluid resistance the particle acceleration. From 
examination the Navier-Stokes equations, they felt that 
the body shape. The group the 
acceleration modulus—was found correlate all the 
values one smooth curve, and Reynolds Number 
dependency was apparent. They felt that their data were 
limited, and planned study the effects decele- 
ration, shape and finite initial velocity over sufficient 
velocity range properly establish the Reynolds Num- 
ber dependency. 

Keim 
Iverson and Balent’s work with circular discs and the data 
both investigations agree very well—a situation which 
unique the reporting acceleration drag 
Keim also measured values for cylinders 
various Jength-to-diameter these showed 
explicit dependency the Re, and 
parameters. Figure the the survey 
have replotted some Keim’s data using and 


has more recently repeated and extended 


clearly show the manner which 
parameter acts displace the curve. 
One noteworthy characteristic which was also indicated 
Lunnon’s data that the effect acceleration the 
drag becomes reduced the Reynolds Num- 
ber, and therefore the turbulent component the wake 
structure, increased. The effect the ratio 


STANDARD 


200 400 600 1000 2000 4000 6000 


Figure 2—Drag coefficients cylinders accelerating from 
rest (Keim). 


the values was very similar that found steady 
motion. Keim’s dimensional analysis derivation the 
acceleration modulus more complete than the one 
given Iverson and Balent that predicts depen- 
dency the rate change acceleration well 
higher order derivatives the velocity which would 
correlated dimensionless parameters the general 
form Keim made attempt, how- 
ever, measure test for the effect the terms 
higher order than the acceleration. 

Working with system similar that studied 
Keim, Bugliarello published lengthy series results 
obtained with 7.6, 8.9 and 15.2-cm. spheres which were 
drawn constant forces through stagnant pool 
water. Again the added mass concept was found 
completely inapplicable, giving uncorrelated added mass 
which ranged from high 500 down 
negative quantities. Bugliarello found that his results were 
best represented the total drag coefficient His 
plots for the 7.6-cm. 
spheres, though considerably scattered, are qualitatively 
similar the right hand section Lunnon’s curves, and 
they show drag coefficients with values high 30, 
compared the steady state value 0.4. The results 
with the 15.2-cm. sphere were considerably more scattered 
and indicated smaller acceleration effect; many cases 
the drag had been reduced almost half their 
steady state values, rather than increased the other 
two cases. Bugliarello considers that this be- 
havior was due the sphere boundary layer being made 
prematurely turbulent surface roughness, and correct 
this again indicates reduced acceleration effect possibly 
brought about increase wake turbulence. 

The authors this survey have replotted some 
Bugliarello’s data again using and co-ordinates 
and the acceleration modulus the varying parameter, 
and the result, shown Figure indicates the same 
diminishing the with increasing 
Reynolds Numbers. The curves roughly show the same 
features the cylinder results Keim, but Bugliarello’s 
case the data are too scattered and limited establish the 
curves quantitatively. 

maintains, and his data seem indicate, that 
values result which can derived from 
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Figure 3—Drag coefficients for accelerating spheres 


(Bugliarello). 


potential flow considerations, but will shown, the 
flow about sphere described potential theory only 
the very beginning acceleration, and the reported 
tendency may have been due the fact that the high 
acceleration data always corresponded the beginning 
the sphere motion. 


The rectilinear but varying acceleration encountered 
particles falling being through stagnant 
fluid represents but one the many possibilities which 
may occur solids-gas systems. For example, 
rest inserted into rapidly moving gas stream would 
experience absolute acceleration until equilibrium 
velocity attained, but considering its behavior relative 
the fluid, instantaneously has finite Reynolds Num- 
ber which diminishes the particle decelerates relative 
the fluid the equilibrium velocity. the other 
hand, particle fired into stagnant fluid from nozzle 
with velocity greater than its terminal velocity would 
experience both relative and absolute deceler ration, and 
extrapolation drag results obtained from the former 
system the latter would not valid, particularly since 
the acceleration drag mechanism probably influenced 
both relative and absolute effects well their 
signs and magnitudes. 

The system consisting particle injected from rest 
into moving gas stream most often occurs the feeding 
zone solids-gas apparatus and its investigation 
ported Ingebo recent series papers deal- 
ing with the acceleration small droplets solid 
spheres injected cloud with negligible entrance 
velocity into air stream moving with velocities 140 
180 feet per second. Under these conditions, the par 
ticles had accelerations large 60,000 value 
far greater than the maximum for particles 
falling into stagnant fluid. The velocities were measured 
fixed stations from downstream from the 
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Figure 4—Drag coefficients spheres entrained 
moving air stream (Ingebo). 


Surprisingly all the data, regardless the acceleration 
magnitude whether they were obtained with solid 
spheres liquid droplets, fell the single line shown 
igure and were correlated the relationship 


Cpa = (5) 


for the entire experimental Reynolds Number 
compared the steady state value 0.55, and this 
trend can shown extend higher Reynolds Num- 
bers, the drag coefficients will fall very small values. 

The particle “drag was not affected the presence 
other particles, since ten-fold increase their concen- 
tration did not bring about any noteworthy alteration 
their velocities. 

results are unique several respects: firstly, 
they show decrease the drag rather than 
the increase due acceleration reported most workers 
for this Reynolds Number region; secondly, they show 
departure from the standard curve the 
Reynolds Number increased, whereas all the other accele- 
ration data discussed far show decreased accele- 
ration effect the Reynolds Number increased; and 
finally all the drag data regardless the magnitude 
the associated accelerations fall single curve rather 
than the individual curves usually obtained for 
ent accelerations, although, this may have resulted from 
the experimental system being one which the Reynolds 
Number and the acceleration could not varied inde- 
pendently and the highest accelerations are here associated 
with the highest Reynolds values obtained. data are 
amount free stream turbulence and 
tainty the accuracy the fluid velocity determina- 
tions; nevertheless, they dramatically underline the neces- 
sity investigating all the possible 
velocity systems individual basis. 

used observe spray droplets entrained and accelerating 
should not differ that the solid spheres. data 
shown Figure indicate clear dependence air 
velocity, unlike those given Ingebo. results are 
surprisingly low, times representing decrease one 
thousand percent from the standard curve values. Fledder- 
man and Hanson felt that this was due evaporation 
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Figure 5—Drag coefficients for spray particles entrained 
moving air stream (Fledderman and Hanson). 


occurring the droplet surface, but this possibility seems 
have been ruled out work, which involved 
both evaporating aporating droplets. The 
particle velocity data may have been error since they 
were not obtained directly, but were derived from 
Froessling’s mass transfer equation relating the drop- 
let velocity the photographically determined decrease 
diameter with time. The application this equation 
sprays droplets has been more less verified 
recent work Manning that the error the 
drag data probably not excessive. 


The large difference between the data Ingebo and 
those Fledderman and Hanson remains explained, 
well the complete contradiction between these results 
which indicate drastic reduction the drag compared 
with the increase caused acceleration the other 
experiments. The fluid turbulence intensity Fledderman 
and Hanson’s test section was found the order 
two percent (based measurements taken when 
particles were present), whereas the turbulence Ingebo’s 
apparatus, while not measured, was probably much 
because the presence 40-mesh screen the 
upstream section. may that the disagreement part 
reflects reaction varying relative intensities and scales 
turbulence, but the information from the experiments 
not sufficient establish any dependency. second 
complication may result from the necessity correlate 
all acceleration drag data function time para- 
meter which related the total time allowed for the 
build-up the boundary layer and wake the particle. 


Rectilinear, Constant Acceleration 

discussed the previous section, the drag measure- 
ments taken using constant driving forces such the pull 
gravity weighted string attached the body, 
well those measured particles moving freely 
through gas stream, are always associated and thereby 
1948 initiated program which removed this compli- 
cation using towing device which could set 
pull the test body different constant acceler ations, and 
the required drag force was instantaneously 
tinuously recorded throughout each run means 
strain gauges. Experiments with circular discs set perpen- 
dicular the motion showed two categories drag 
behavior. the first the increase drag remained 
roughly proportional the acceleration but 
the second, the drag was initially much larger but 
eventually dropped the values found the first case. 
Insufficient data were reported indicate any explicit 


VERTICAL CYLINDER 


Ud/U, 


Figure 6—Drag coefficients cylinders and spheres sus- 
pended oscillating fluid (Crooke). 


Reynolds Number time dependency. For circular 
cylinders, however, was actually less than for 
10%, but surpassed the region 
the increase due the acceleration increasing 
cylinder diameter decreased, and the region about 105 
Luneau was unable obtain reproducible results. Laird and 
Johnson 1956 also reported long-range program 
study the acceleration effects the drag coefficient, 
and their preliminary results with cylinders 
length towed constant acceleration through water did 
not show any Reynolds Number correlation, and only 
showed slight correlation with the 
meter. The data are very scattered, but they generally 
show that for the Reynolds Number range 1.5 
1.25 105, went from 0.6 1.7 (Ud,/U,?) went 
from —0.04 0.18 compared with the steady state value 
0.85. These data are among the first result from 
studies the effects deceleration, and the reduction 
the drag coefficient this case particularly note- 
worthy. Data from constant acceleration experiments are 
just beginning accumulate that even qualitative con- 
clusions cannot yet made, other than that the results 
obtained with varying acceleration differ from those 
which the acceleration held constant. 


Acceleration Oscillatory Motion 

The post-war period has given rise the renewed 
investigation the acceleration effects particles oscil- 
lating fluid. These experiments are divided into two 
groups: the first, the body held fixed while the 
ambient fluid given specific oscillation—in some cases 
resembling the motion waves fixed pier; the 
second, the test body moves and fro usually attach- 
ment system springs. The important consideration 
that both systems the continual reversal the rela- 
tive flow direction does not allow the flow patterns 
around the particle develop they would 
linear motion, and the particle continually enters flow 
regime which has been disturbed the previous cycle; 
nevertheless, the data obtained light 
acceleration phenomena general, and will 
larly useful the consideration the behavior small 
particles imbedded large scale turbulent flows. They 
are also important explaining results being obtained 
with solids-gas systems subjected sonic energy. 
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(i) Fixed Bodies Oscillating Flows 


The forces bodies immersed and held fixed 
oscillating fluid have been conventionally correlated 
modified version Equation (1), hich has the form: 


where the positive numerical value and the 
coefficient generally taken being equal 
k,), the addition the added mass coefficient 
resulting from the fact that this system the entire fluid 
—rather than the body—is accelerating, and co- 
efficient which not confused with either the 
steady state the total acceleration drag coefficient 
This equation gives the erroneous impression 
being theoretically well founded, but actually based 
intuitive extrapolation from inviscid fluid behavior 
which assumes that the drag forces during unsteady 
motion finite Reynolds Numbers can parceled off 
into two discrete packages which are independently related 
the acceleration and velocity, respectively. also gives 
rise the possibly unsound practice solving for two 
experimental force coefficients from only one experi- 
mental force determination, thereby complicating the 
interpretation the data. 

Very often the values and are obtained 
considering the oscillation sinusoidal, that the 
term zero when the value drops 
zero the phase angles 90° and 270°, and similarly the 
term eliminated and 180° when the value 
becomes zero—needless say, the values and 
obtained and 90° respectively not generally corre- 
spond with those indicated 270° and 360° they would 
the concept was valid. other instances, the value 
employed the calculation the average 
obtained during the whole series trials, and third 
practice uses the added mass coefficient predicted 
inviscid flow theory. 

series papers, Morison and O’Brien 50, 
conjunction with others, reported results with fixed 
cylinders and spheres, and Equation (6) when applied 
the results for cylinders showed variations from 
1.4 1.6. was found fall with increasing Reynolds 
Numbers from about 1,000 close the 
steady state value 10,000, and there was consider- 
able scatter the data, varying much 700 
percent constant Reynolds Number. Crooke 
analyzed Morison’s data, plotting the total drag coefficient 
The results for spheres and cylinders immersed oscil- 
lating fluids are shown Figure and the correlation 
surprisingly good, showing considerable improvement 
over that given Morison. Crooke found clear 
Reynolds Number trends for the range covered but 
the lower values the acceleration modulus rough quali- 
tative tendencies which are those found with 
rectilinear motion are indicated. about the same time, 
McNown presented paper dealing ith preliminary 
investigations the forces flat plates and lenticular 
cylinders fixed fluctuating flow, and proposed 
Equation (6) that the added mass coefficient ac- 
knowledged being function time, giving the 
relationship: 


McNown reasoned that since and both depend 
the particle shape there must unique relationship 
between them, and single curve with and 
coordinates was theoretically calculated from extension 
the Kirchhoff approach discussed Part this 
series, and was incorporated with artificial model 
the flow which was necessary ignore the com- 
plexities the wake mechanics. Assuming the curve 
valid, then tried out different variations with 
time until one was found which predicted the variations 
force observed. The degree with which the 
relationship corresponds reality difficult establish, 
particularly since the remainder the analysis into which 
inserted may not for real fluids. 

Wiegel, Beibe and Moon measured the forces 
cylinders submerged actual ocean waves whose velo- 
cities were experimentally recorded. spite Crooke’s 
success with the overall drag coefficient the authors here 
again analyzed their results the use Equation (6) 
and obtained very large number points, which were 
unfortunately accompanied very wide scatter 
coefficients showi ing decrease well increase 
with acceleration. The authors felt that the results 
differ from those obtained laboratory trials because 
the variable turbulence present the ocean waves. 
conclusion, they re-analyzed small portion their 
data using the total drag coefficient and here obtained 
moderate success and the results fell between the curves 
given Crooke and Keim, but they unfortunately did 
not have time re- analyze the bulk ‘of their data this 
way. 

Borgman has developed equations vectorial and 
matrix notation which would the drag force data 
already obtained applied cylinders which were 
oriented any arbitrary manner with respect the 
oscillation direction, but the results must checked 
experiment, since the method implies alteration the 
fluid dynamic behavior due the orientation. 

Very recently, Keulegan and Carpenter what 
perhaps the finest work its kind, measured the forces 
cylinders and plates fixed oscillating fluid. They 
modified Equation (5) the form: 


where the instantaneous drag force and the 
cients are here the averaged values taken over the entire 
makes for the fact that and are averaged rather 
than instantaneous values. very important contribution 
made this paper the introduction the fundamentally 
meaningful reduced time parameter U,, being 
the maximum fluid velocity the cycle and its period, 
and this parameter proportional the number 
diameters travelled the fluid one direction before 
the flow reverses, that reflects how much time 
given the boundary layer and wake develop. This 
information may the oscillating flow problem 
the frequency the turbulent fluctuations the case 
small particles imbedded very large scale turbulent 
fields. 

The applicability the time parameter illustrated 
the remarkable degree which correlated the 
resistance coefficients. Smooth curves and vs. 
(U,,T/d) were obtained for both cylinders 
with the results for the various trials all falling close 
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Figure 7—Drag coefficients flat plates and cylinders 
fixed oscillating fluid (Keulegan). 


single curve for each shape. Reasonable explanations 
were given for the different curves obtained for 
cylinders and for plates, but this could not, however, 
done for the added mass coefficients, which should have 
remained constant. Although the authors 
main emphasis the Equation (8) type analysis, plots 
were also given function the time para- 
meter, and the correlation was excellent. shown 
Figure the data for cylinders and plates both follow 
similar curves which the coefficients decrease rapidly 
first and then more gradually the steady state value 
the time parameter increased. Surprisingly absent 
Reynolds Number effect, but wider experimental range 
including different bluff body shapes will now required 
confirm the uniqueness the group 
correlating parameter for oscillating motion. discussion 
this parameter continued the following section 
the fundamental aspects the acceleration problem. 


Oscillating Bodies Fluids 

The post-war investigations the force system 
body oscillating fluid rest have not been actively 
pursued the category just discussed. Carstens using 
table-tennis submerged various viscous fluids and 
attached oscillating weights, felt that 
could best describe the forces the sphere expres- 
sion which equivalent Equation (1), except that 
the (U,?) term was replaced viscous damping force 
which was assumed proportional rather than 
U,?, resulting the expression: 


where “viscous damping” coefficient. Following 
analysis, Carstens led believe that both 
and the ratio the viscous drag predicted 


Stokes Law (the ratio being given the symbol N,) 
where the angular velocity the which 
assumed sinusoidal and 2E/d where the 
maximum amplitude the oscillation. was found 
decrease and the added mass coefficient was found 
increase with increasing values N,, but more experi- 
mental data covering broader experimental range are 
required before any generalizations can made. 


Stelson felt that the scattered added mass values 
reported the literature were more the result inap- 
propriate measurement rather than breakdown the 
added mass concept. Stelson’s experimental technique con- 
sisted attaching the submerged test bodies rigid 
bar which was then vibrated the natural fundamental 
frequency the combined system. The system was first 
vibrated air using varying masses for test bodies 
that relationship between the body mass 
frequency the vibration was found. The added mass 
given body was obtained vibrating when sub- 
merged water and then noting the mass corresponding 
the new vibration frequency. Stelson ascribed the 
apparent increase mass being entirely due the 
added mass the body. The amplitudes the oscilla- 
tions these experiments are very much smaller than 
are usually obtained experiments with oscillating bodies, 
that the fluid system about the body does not have time 
depart from the potential flow pattern, and for this 
case alone the classical hydrodynamic treatment may pro- 
vide good approximation. 

Stelson obtained results for variety shapes well 
spheres, and for the latter, was found equal 
0.51, which almost coincides with the theoretical predic- 
tion. For cubes, was found equal 0.67 for 
both the broadside-on and edge-on motion, and although 
there was theoretical solution available for the cube, 
Stelson points out that and had 
shown that the theory predicts that the added masses 
both positions should agree with each other. The values 
for rectangular with varying length-to-width ratios 
and for circular cylinders with 
diameter ratios fell smooth curves extrapolating 
values predicted potential flow theory. The results 
obtained, together with the analysis employed, seem 
valid for the one extreme the oscillation problem 
which the amplitude negligible. The experiments 
previously cited point out, however, that with increasing 
amplitude the problem rapidly enters into regime 
which the empirical and analytical difficulties reflect the 
complexity the attendant flow characteristics. 


Fundamental Investigations Accelerating 
Phenomena 

The following discussion divided into two sections: 
the first deals with attempts elicit information con- 
cerning the changing flow fields and drag coefficients 
unsteady motign approximate solution the pertinent 
Navier Stokes Equations, and the second discusses the 
empirical investigation the effect non-steady motion 
the boundary layer and wake properties. 


Theoretical Investigations 
(i) Reynolds Numbers 

the system confined one which the accele- 
rations are very moderate, and the particle Reynolds 
always less than one, becomes possible 
overcome the difficulties introduced the non-linear 
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Figure 8—Drag coefficients for cylinder impulsively 
accelerated from rest (Payne). 


terms the Navier Stokes equation, and calculations 
Basset Boussinesq and Oseen show that the 
resultant drag force particle starting from rest and 
accelerating stagnant fluid given by: 


The first right hand term identical the added mass 
resistance predicted classical hydrodynamics and the 
second term coincides with the steady -state drag given 
Stokes Law. The third term shows that even the 
idealized Stokes Law region the drag force will depend 
the previous history the particle motion. Tchen 
has generalized Equation (10) for systems where both 
the particle and fluid are undergoing acceleration, but the 
forces will become complex functions the Reynolds 
Number well the acceleration parameter (U,d/U,?). 

More recently, Pearcey and Hill have dealt with 
the case spherical particle which instantaneously 
Their solution employs the Oseen approximation 
the inertia terms (see Part that limited 
Reynolds Numbers less than unity. For small values 
the time parameter the drag force 
given by: 


and for large 


The solutions indicate that the drag force impulsive 
acceleration made three components: the first 
identical with the steady state Stokes Law drag, the 
second infinite the start the motion and decreases 
and the third component acts only the instant 
that the impulsive motion Equations (11) and 
(12) suggest that the region extremely low Reynolds 
Numbers the effect the non-steady motion the drag 
coefficient will decrease with decreasing fluid densities, 


and specific solutions for the trajectories small spheres 
air should show only small effect due the accele- 
ration compared the large effects noted experimentally 
air higher Reynolds Numbers. 


(ii) Reynolds Numbers 


There does not seem any overall unsteady-state 
drag solution for the post-wake-forming Reynolds Num- 
bers for the sphere system, that, before, recourse 
made work done for the simpler case the cylinder. 
Payne using finite difference numerical approach 
similar the methods employed Kawaguti and 
Thom for the steady state cylinder case (see Part 
obtained specific solutions for impulsive motion occur- 
ring Reynolds Numbers and 100. The solutions 
are plotted Figure showing the drag coefficient 
unique function the Reynolds Number and the 
time parameter (U,6/d), the latter being equal the 
number diameters traversed the cylinder since the 
initiation the motion. The drag drops rapidly first 
and then gradually the steady state value surpris- 
ingly short period time, requiring only distance 
about three diameters come within the value 
given the standard curve. the beginning the 
motion the ambient boundary layer very thin and the 
intervals employed the numerical computation are too 
large obtain accurate values For this region 
Payne reverts boundary layer solution offered 
Goldstein and Rosenhead which predicts that for 
tending zero: 


the skin friction accounting for half the total. This 
boundary layer solution and the others, which will soon 
discussed, are not reliable after separation occurs (for 
the cylinder this happens when approximately 0.2 diameter 
has been traversed) that Payne has joined the valid 
regions both solutions with dashed line. note- 
worthy that the solutions show dependency the 
fluid density, beyond that contained the Reynolds 
Number—in contrast with the equations previously cited 
for Reynolds Numbers less than unity—and the total 
drag coefficient was easily represented single smooth 
curve. Payne expects extend his treatment higher 
Reynolds Numbers the use improved electronic 
compared with empirical measurements already available 
for Reynolds Numbers the vicinity 600. 

attempt has been made Schwabe obtain 
analytical solution for the pressure fields about impul- 
sively -accelerated higher Numbers, 
using potential flow theory. The results were compared 
with experiment, and the case with steady motion 
this approach was moderately successful for the frontal 
region the wake, but the wake region had 
replaced artificial model allow the calculation 
the non-wake region. the time increased, the dis- 
crepancy between the calculated and observed flow pat- 
terns became larger and the region applicability 
potential flow analysis was shown restricted the 
beginning the motion. 


(iii) Boundary Layer Solutions 

The boundary layer solutions available for impulsive 
motion are valuable two respects, that they describe 
the fluid motion the very beginning the event, and 
they predict the time interval between the initiation 
motion and the separation the boundary 
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Blasius launched the mathematical attack the 
problem the boundary layer build-up cylinder 
impulsively accelerated from rest his approximate 
solution indicated that separation will occur after the 
cylinder had moved the short distance 0.15 diameter. 
describes the initial few moments being transfor- 
mation from potential state flow with boundary 
layer zero thickness developing initially through fric- 
tional processes and subsequently through the convective 
terms the equations motion. The boundary layer 
soon thickens and separates, and important realize 
that the boundary solutions are longer valid past 
this point. For the case constant acceleration from 
rest, Blasius found that separation occurred when only 
0.6 diameter had been traversed that generally the 
period application boundary layer solutions 
particularly short one. The wake build-up shown 
occur almost the start the movement interval 
which would not detected 
Improvements the Blasius solutions have been made 
and Rosenhead but the corrections the 
time intervals for separation are 

Boltze provided approximate solution the 
boundary layer equations for impulsively -accelerated 
sphere, and the results were similar those obtained with 
cylinders with separation occurring after 0.2 diameter had 
been traversed. More recent works (74 75, 76, 77, 78, 79, 80) 
have provided refined solutions which include treatments 
with arbitrary time-dependent velocities with particular 
emphasis being placed the flow near the front stagna- 
tion point, and they emphasize the rapidity the bound- 
ary layer build-up the point where separation first 
occurs. 


(iv) Solutions for Oscillating 


The analytical treatment body oscillating 
fluid simplified the motion restricted instantan- 
eous Reynolds Numbers which are much smaller than 
one, and for this case Lamb has derived formula 
giving the value the added mass coefficient as: 


where the frequency oscillation. For finite Rey- 
nolds Numbers the flow system becomes far more com- 
plex and the treatment has been thus far usually restricted 
oscillations with amplitudes which are smaller than the 
characteristic diameter the immersed body 
systems which vortex shedding does not occur. The 
boundary layer solutions for bodies fixed oscil 
fluid are closely associated with the experimental observa- 
tions the authors concerned that these will taken 
the experimental section. Solutions for the boundary 
layers cylinders oscillating on-coming fluid have 
been the frontal region the bodies and have 
been given Glauert for oscillations transverse 
the stream and Lighthill for oscillations parallel 
the stream. 


Experimental Investigations 

One the most valuable insights into the processes 
occurring during impulsive acceleration has resulted from 
investigation Schwabe Using minute particles 
trace and record the velocity and direction movement 
the components the flow, succeeded obtaining 
series-of photographs showing the progressive evolution 
the ambient flow structure from the very beginning 
the motion. These tracer velocity measurements were 
sufficiently accurate allow the calculation the entire 
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Figure 9—Drag coefficients and wake dimensions for 
cylinder impulsively accelerated from rest (Schwabe). 


pressure distribution occurring given instant, which 
not possible with pressure probe because the finite 
time required for even single measurement. Integration 
the pressure field yielded the acceleration drag 
cients and these, together with the dimensions the 
various wake components for constant Reynolds Number 
580, are shown Figure function the time 
parameter which represents the number 
cylinder diameters traversed from the start the motion. 
would seem that the drag coefficient rises rapidly from 
zero first and then more gradually maximum which 
almost twice the steady state value. may inferred 
that the eventual decrease Cp, the steady state 
probably very gradual one, although this was not ascer- 
tained Schwabe, since the experimental system did not 
allow measurements beyond There are 
several noteworthy discrepancies between this result and 
the analytical solutions cited previously. Although the 
latter predict infinite drag descending rapidly 
the flow develops, Schwabe’s results indicate the reverse. 
This may partly due Schwabe’s result giving only 
the form drag, omitting the skin friction component. The 
form drag, according the boundary layer solution, 
should also infinite but paralleling the case 
Payne’s calculations, the computational intervals used 
the evaluations the pressure are too large detect 
the gradients occurring the start the process when 
the boundary layer extremely thin. More important— 
but more difficult the fact that Payne’s cal- 
culations suggest the period increased 
cient occurs interval which not greater than 
(U,6/d) 0.5, for and 100, the interval 
well the drag increase due acceleration tending 
decrease with increasing Reynolds Numbers. Schwabe’s 
experimental force determinations— which can only 
conservative due the omission the skin friction com- 
ponent—show much longer period abnormal drag 
behavior well larger increase the drag coefficient, 
this being more compatible with the results obtained 
the macroscopic investigations. more direct comparison 
will available when Payne obtains results higher 
Reynolds Numbers. 

Recalling that Ingebo’s drag determinations were car- 
ried out under conditions which initially approximate 
impulsive acceleration, there should some agreement 
with the fundamental investigations Payne and 
Schwabe. Particularly difficult explain the decrease 
drag shown Figure Although Schwabe does record 
decreased drag the initiation the motion (this trend 
being subject revision when the boundary layer con- 
tribution taken into account more fully) the 
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limited the time taken for the particle travel one 
diameter (at the most), after which increased drag forces 
are shown. 

The picture further complicated reanalysis 
Ingebo’s work given Mercier according which 
the coefficients and Ingebo’s correlating equation: 


Cpa = ky/Reks (15) 


are not equal and 0.84 respectively, had been 
reported Ingebo, but vary widely, having values 
which range from 0.5 Since Ingebo’s apparatus 
there were four observation stations placed downstream 
from the injection point (designated 0,1,2,3), Mercier 
integrated the equation motion which incorporated 
Equation (15) and solved for and using either 
the consecutive intervals 0-1, 1-2, 1-2, 2-3. This pro- 
cedure somewhat inconsistent since the one hand 
assumes that the coefficients remain constant over the 
interval integration, whereas was shown that 
given experiment would sometimes have the value 
0.5 when obtained from the interval 0-2, and the value 
for the interval 1-3, which leads the suggestion that 
may vary continuously throughout the interval 0-3, 
possibly having values less than 0.5 and greater than 
Nevertheless, Mercier suggests that the accelerating par- 
ticles are initially dynamic regime which 0.5 
and they then enter into region which 2.0. The 
transition between the two regimes supposedly occurs 
abruptly points whose loci describe the curve: 


Furthermore, Mercier finds that the drag coefficients 
the beginning the motion fall the curve: 


and explains this being result the flow structure 
around the particle when being impulsively accele- 
rated. This assertion contradicted the available 
fundamental studies since Equation (17) predicts rapidly 
increasing lowering drag coefficients the Reynolds 
Number increases above 40. Although Mercier’s analysis 
seems fit results may that the results 
are subject systematic errors (as was suggested 
Ingebo himself) and complicating phenomena other 
than acceleration, that may too early attempt 
tionships without the further accumulation experi- 
mental data. 

series fundamental studies the wake system 
flat plate accelerating constant rate from rest was 
started Batailler 1956 with the aid 
elaborate photographic technique which allowed fine-scale 
observations taken the velocity and acceleration 
the wake velocity field function time. This was 
accomplished inserting minute tracer particles into the 
flow and then illuminating them with closely spaced 
series xenon lamp flashes known time duration. 
the start the motion, the flow this investigation was 
also found potential but the velocity increased 
there was growth stagnant fluid the edges the 
plate and eventually periodic discharges eddies de- 
veloped. second study, stronger accelerations were 
employed, and here there appeared strange chain 
discharge from the plate edges small stagnant pools 
fluid which eventually fed into the large standing eddy 


the plate rear. The drag coefficient was found 
more than twice the steady state value. 

Using dimensional analysis, Batailler next attempted 
find parameter which could used characterize 
the wake dimensions unsteady system, and this led 
the development the dimensionless time group: 


The changing wake length taken the distance 
from the rear face the plate the rear separated 
stagnation point) was very well for the range 
accelerations and plate diameters the 
relationship: 


Lw ‘d = 0.5 Os. 


with surprisingly little deviation the data this 
straight line relationship. These investigations are yet 
preliminary but the time-parameter should prove 
very useful correlating the accumulating data wake 
dimensions and drag forces with constant acceleration 
rates. interesting note that some ways the 
finite acceleration ana!ogue Schwabe’s impulsive accele- 
ration time-parameter, (U,6/d) which equal the 
distance, body diameters, traversed since the start 
the motion. the simplified case constant acceleration, 
the relationship distance traversed 1/2 can 
substituted into the time group given Equation (18), 
the latter then found proportional the square 
root the distance, body diameters, traversed since the 
start the motion. Similar substitution shows that the 
acceleration constant, the acceleration parameter 
inversely proportional the number body 
diameters traversed, that the constant acceleration 
case alone correlation and becomes 
synonymous, although different curves will obtained. 


exacting study the vortex system formed during 
the oscillation body fluid has been made 
Andrade the course investigation 
into the use solid particles means measuring 
the modes the sonic vibration gases the familiar 
Kundt’s tube method. Both cylinders with 
diameters low microns were individually placed 
into column gas which was vibrated various fre- 
quencies and amplitudes. The flow systems surroundin 
the suspended spheres were made visible the plane illu- 
mination minute smoke particles, and when the Reynolds 
Number the oscillating gas flow entrained (calculated 
using the maximum relative velocity U,,) exceeded 
critical value there was clearly revealed not one but two 
vortex rings opposite ends the plane vibration. 
Oddly enough, the flow the vortices did not oscillate 
but was steady, and furthermore circulated direc- 
tion opposite that which occurs the steady state case. 
Andrade found that for the larger spheres employed (0.2 
cm. 0.6 cm.), the Reynolds Number for wake 
formation, was function the ratio the oscilla- 
tion amplitude the particle diameter regardless 
the frequencies (see Figure 10). For the small spheres 
(0.025 cm. 0.042 cm.) there was dependency 
the (E/d) ratio all, even though was varied over 
large range, and Re, had constant value which 
surprisingly close the value obtained for steady flow. 

Assuming that the results are not affected support 
interference, the different behavior noted with small par- 
ticles compared with large ones raises very important 
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point regarding the application data obtained from 
experiments with the latter the small particles occurring 
solids-gas systems. Andrade offers possible explana- 
tion for the difference referring mathematical 
solution for the boundary layer thickness oscillatory 
motion given Lord Rayleigh noting that the 
ratio the boundary layer thickness calculated from 
this solution the partic cle diameter much larger for 
the small spheres than for the large ones. This sug- 
gests that unsteady phenomena observed with large 
bodies must small particle experimenta- 
tion before can made. 

which was similar Andrade’s, obtained photographs 
which shwwed the same vortex formation, and suc- 
ceeded providing boundary layer solution 
successive approximations which was 
agreement with Andrade’s results, but also predicted 
‘an inner vortex system the immediate vicinity the 
cylinder which moved the opposite direction the 
outer one. Its presence had not been observed either 
Andrade Schlichting himself. about the same 
time, Carriere 92) studying similar system found that 
the vortex which observed corresponded the ones 
found Andrade and Schlichting, but the direction 
movement was reversed. 


More recently there has been flurry activity 
workers acoustics seeking more accurate description 

“acoustic streaming”. Skogem and West continued 
Andrade’s experiments with cylinders and detected the 
presence the inner vortex Schlichting. 
Westervelt concluded that the direction flow 
the outer vortex depended the amplitude-to-body 
diameter ratio (E/d), while Andres experimenting 
with cylinders found that low values the Reynolds 
Numbers (based maximum relative velocity only the 
inner vortex system appeared, but the outer flow devel- 
oped when the (E/d) ratio was increased. Andres 
proceeded solve the boundary layer equations for 
large and then small Reynolds Numbers (using the Oseen 
the latter) and the observed flow 
patterns were indicated the solutions for low values 
the velocity amplitude. analysis suggested that the 
type flow occurring was characterized the group 
(Re The solution using the Oseen approximation 
may not valid since esterveld (38) has developed 
proof showing that the use Oseen’s approximation 
derive the flow field the vicinity oscillating object 
is, the words the author, 

Experimental work with cylinders has been continued 
and with spheres Lane Skavlem and 
improved the solution the Stokes equa- 
tions and obtained very good agreement with Holtsmark’s 
observations. Lane has modified the method solution 
deal with spheres and showed that there was good 
similarity between the properties the two flows. his 
experiments with spheres was careful check for pos- 
sible suspension interference and was able avoid this 
limiting his diameter range minimum 3/16-in. Al- 
though there was fair agreement between experiment and 
his theoretically predicted boundary layer thickness, closer 
agreement streamlines—especially the lower oscilla- 
tion frequencies has since been obtained Thrasher and 
Westervelt who have altered the theoretical treat- 
ment. 


LARGE 
SPHERES 


SMALL 
12 SPHERES 


0.00! 0.01 


Figure 10—Minimum Reynolds Number for vortex for- 
mation for spheres fixed oscillating air column 


(Andrade). 


The above work vividly points out the differences 
between the wake structures occurring oscillating and 
steady motion, and similar differences which must influ- 
ence the drag coefficient most likely occur the pre- 
wake-forming Reynolds Numbers. some the experi- 
ments there was steady motion superimposed over the 
oscillating fluid motion. This caused the upstream ring 
diminish while the downstream vortex was distorted with 
increasing mean fluid velocity. Although the oscillation 
system gradually obliterated with increasing 
mean motion still exerts important influence, 
turn distorting the steady state wake. Measurements 
the drag forces all these experimental studies, would 
have been very valuable but they could not taken 
without difficult modifications the equipment. 

Very little known about the effect the oscillations 
the lower critical Reynolds Number which the 
wake shedding first begins. Carstens seems have 
uncovered rather profound influence rough check 
made with hot-wire anemometer during his oscillating 
table-tennis ball experiments which showed that vortex 
separation had occurred even when the maximum Rey- 
nolds Number during the oscillation was high 3600. 
Schwabe’s work indicates that this effect must related 
the period the oscillation, the critical Reynolds 
Number increasing with the frequency. Keulegan and 
Carpenter and McNown and Keulegan suggest 
vortex development oscillating systems should cor- 
related the ratio the time flow one direction 
the time taken form vortex the steady 
state case, the latter being obtained from the Strouhal 
Number (see Part and they show that the para- 
meter successfully correlated their data 
—is related the (T/T,) ratio through the relation: 


(T/T,) 


repeated vortex shedding does occur unsteady 
flow, then the corresponding Strouhal Number most 
likely varies throughout the cycle and its relationship 
the steady state would not simple one. McNown 
and Keulegan found that the fluid distortion the 
previous cycle inhibited the vortex formation and they 
underline the need for extensive investigation the 
wake properties unsteady systems. 
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Conclusions 


Although the available experimental and theoretical 
information concerning the unsteady motion particle 
through fluid more deficient and cruder condition 
than that which available for the corresponding steady 


the following itemized generalizations are 


suggested: 


(a) 


(b) 


(e) 


(f) 


(g) 


potential flow theory can used describe events 
only the beginning rectilinear accelerations 
oscillatory motions involving very small 
amplitudes, and the added mass concept loses theo- 
retical significance and practical utility applied be- 
yond this region; 

non-steady motion can cause very appreciable de- 
partures from the steady-state drag even for particu- 
late movement gases where 
density ratio extremely low; 


acceleration phenomena are correlated reduced 
time parameters related the number particle 
diameters traversed since the initiation the motion 
and are also dependent the Reynolds Number 
describing the flow; 

the use the total drag coefficient gives better 
results and has more theoretical justification than the 
summation coefficients often employed; 

constant and varying acceleration rate experiments 
not give similar drag coefficients and therefore the 
rate change the acceleration must taken into 
account; 

increase wake turbulence caused Reynolds 
Number increases surface roughening seems 
decrease the effect the acceleration the drag; 
and 

the large decreases drag noted some accelera- 
tion experiments can not explained terms 
present fundamental information 
increases drag only and suggested that these 
data are influenced phenomena other than accele- 
ration. 


From what being said the workers participating 
the recent renewal activity this field, there will 
many energetic inquiries made with view pro- 
viding reliable quantitative drag values well further 
insight into the fluid namics unsteady motion. 
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Nomenclature 


Any consistent set units may employed. 


Those listed 


are merely illustrative. 


wake length, ft. 

steady state drag coefficient, dimensionless. 

total drag coefficient accelerated motion, dimension 
less, Equation (4). 

drag coefficient, dimensionless, Equation (6). 

inertia coefficient, dimensionless, Equation (6). 

body diameter, ft. 

oscillation amplitude, ft. 

oscillation frequency, cycles/sec. 

body force due gravity, lb. force. 

local acceleration due gravity, 

distance between vortex centres, ft. 

added mass coefficient, dimensionless, Equation (1). 

coefficient, force Equation (2). 

coefficient, (Ib. 2)/ft., Equation (2). 

viscous damping coefficient, (Ib. 

tion (9). 


Equa- 


= 


coefficient, dimensionless, Equation (15). 

coefficient, dimensionless, Equation (15). 

cylinder length, ft. 

characteristic particle length, ft. 

wake length flat plate, ft. 

axial distance between vortex centre and rear stagnation 
point, ft. 

particle mass, 


added mass, Equation (1). 


oscillation parameter pw/d, 

viscous damping ratio dimensionless. 

exponent, dimensionless. 

steady state fluid drag force, lb. force. 

acceleration fluid drag force, force. 

supplementary drag force, force. 

particle Reynolds Number, dimensionless. 

minimum Reynolds Number for wake formation, dimen- 
sionless. 

Strouhal Number, dimensionless. 

period oscillation, sec. 

time parameter 

maximum relative velocity oscillating cycle, ft./sec. 

fluid velocity, ft./sec. 

body velocity relative fluid, ft./sec 

positive numerical value U,, ft./sec 

particle volume, 

linear distance, ft. 

= time, sec. 

reduced time parameter dimensionless. 

reduced time parameter ‘d, dimensionless. 

fluid density, 

angular velocity, radians/sec. 
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FORMULA” INSULATIONS ARE SUBSTITUTE 
FOR EXPERIENCED APPLICATION ENGINEERING! 


Experience has proven that for correct, economi- 
cal insulation there still only 
approach: experienced application engineering. 

Canadian Westinghouse continues offer the 
widest range motor insulations meet 
every application need. These insulations have 
been applied all classes Westinghouse 
motors, and their properties tested under the 


most varied conditions industry. Among 


advanced Westinghouse motor insulations are: 


ENCAPSULATED INSULATION 


armament system for mush-wound motors. The 


a new 


stator completely encapsulated epoxy. The 
method used apply the epoxy gives maximum 
cooling and protection the coils. Encapsulated 
insulation for conditions subject 
excess humidity industrial concentration 


where abrasives other dusts collect motor 
windings. 


THERMALASTIC INSULATION with the 


highest insulating qualities known for large 


motor stator THERMALASTIC INSU- 
LATION has high thermal stability and effective 
heat dissipation. completely free voids, 
and retains its elastic strength. Resists moisture, 
carbon black, dirt, oils, acids and alkalies. 
mica-resin insulation, and the resins can 
used blocks.’’ This enables 
Westinghouse scientifically design the parti- 
cular THERMALASTIC INSULATION that 


most suitable for every class service. 


Insulation key factor the cost industrial 
drives. standard service 
severe conditions, are both 
ineffective and costly. Maximum economy 
possible only when the insulation selected 


suitable every way the conditions service. 


Westinghouse the original pioneer the 
development motors with more experience 


motor insulation than any other manufacturer. 
Today, always, Westinghouse motors are 


engineered, built—and insulated—to meet your 


applications best! 


you CAN Westi nghouse 
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packaging and electric frying pans, Everyone’s using aluminum, 


Why not our 


Foil packaging that seals freshness and flavour 
appliances that are light, strong and 
convenient. Aluminum everywhere. 
wonder. other metal offers such remarkable 
combination qualities. And its uses continue 
increase all the time with the development new 
alloys, improved fabricating and welding techniques 


—and growing consumer 
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“KNOW HOW” CAN HELP YOU YOUR BUSINESS. 
ALCAN are the people see about everything con- 
cerning aluminum. They are leaders its develop- 
ment and set its standards quality. ALCAN has 
over fifty years experience aluminum and the 
major source Canada for aluminum sheet, wire, 
rod, bar, tubing, foil, extrusions, castings and 
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